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We derive an equilibrium statistical theory for the macro- 
scopic description of a ferromagnetic material at positive finite 
temperatures. Our formulation describes the most-probable 
equilibrium macrostates that yield a coherent deterministic 
large-scale picture varying at the size of the domain, as well 
as it captures the effect of random spin fluctuations caused by 
the thermal noise. We discuss connections of the proposed for- 
mulation to the Landau-Lifschitz theory and to the studies of 
domain formation based on Monte Carlo lattice simulations. 
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I. INTRODUCTION AND SUMMARY OF 
RESULTS 

Based on a lattice model recently used for micromag- 
netic simulations at positive finite temperatures (see |l|]) 
we derive a statistical equilibrium theory that describes 
large-scale coherent structures (magnetic domains, do- 
main walls) in the presence of thermal fluctuations. The 
resulting description of equilibrium states depends only 
on macroscopic quantities: magnetization m, induced 
magnetic field h, external magnetic field as well as 
structural properties of the ferromagnetic material (crys- 
talline anisotropy). Furthermore, our formulation cap- 
tures thermally induced spatial, as well as spin random 
fluctuations on the lattice and their effect in the macro- 
scopic model. The Landau-Lifschitz theory (see [§,||) has 
been successfully applied to modeling of magnetic mi- 
crostructures in ferromagnetic materials with crystalline 
anisotropy. However, the assumptions under which the 
theory is derived rule out systematic description of ther- 
mal noise and the theory is sometimes viewed as a zero- 
temperature approximation. Recent advances in applica- 
tions of ultra-thin magnetic films pointed out the impor- 
tance of thermal fluctuations for studying such effects as 
spin reversal, nucleation, metastability or hysteresis (see, 

e.g., i, i). 

The approach proposed in this communication pro- 
vides systematic treatment of the micromagnetic lattice 
model at a (positive) finite temperature It allows for the 
derivation of, (i) the probability distribution ( 



of 

spins at all points of the magnetic specimen by incorpo- 
rating the entropy, derived fro m the microscopic model, 
into the continuum functional (1.12); (ii) a macroscopic 
Ginzburg-Landau type free energy (1.9) for the average 



magnetization; (iii) the description of random fluctua- 
tions in configuration space given by (1.8) and the derived 
large deviation principle ( |2.5| ) respectively. In the case 
of a large specimen where one observes well-developed 
magnetic microstructure with many magnetic domains, 
our formulation also provides an algorithm to compute 
spatially averaged magnetization fields. In a follow-up 
publication we study stochastic dynamics and domain 
formation and evolution in the context of the statistical 
framework developed here. 

The physics on the regular lattice £jv C i?'' (d = 2, 3) 
with TV sites is defined in the standard way by the 
means of an interaction potential between spins cr(a;), 
a{x') with values on the unit sphere S"^ at two differ- 
ent sites X, x' G Cn- Since we do not pursue an ab 
initio derivation, we assume that the interaction poten- 
tial U{x~ x', cr{x), o'(x')) consists of the following terms: 
U = Ue + Ua + Ud + Uh with particular contributions re- 
flecting different types of interactions between magnetic 
moments. The underlying crystallographic lattice struc- 
ture is incorporated through the interaction energy den- 
sity 4" rather than by using different lattice geometries. 
The exchange energy is given by 

= -| E E ^((^ - ^')/N'/')<j\x)a^ix) , (1.1) 

x,x' i,j 

with the local mean-fleld interaction described by the 
positive function J. This term is often approximated 
by the nearest-neighbor interaction only. The anisotropy 
energy related to the crystalline structure of the material 
is defined by the energy density ^E*, 



Ua^Kj2^{a{x)). 



(1.2) 



The non-local, long-range interaction between different 
magnetic moments (spins) is described by the dipole- 
dipole interaction 
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a'{x)cr^x'). (1.3) 



With a slight abuse of notation we write a{x)a-{x') = 
j cr^{x)a^ {x') for the scalar product of two unit vec- 
tors at sites x and x' and points on S*^ are represented 
by unit vectors with the components cr\ i — 1,2,3. We 
define the interaction Hamiltonian of the system 
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HN{<y) = -^ ^ U{:x - x ,a{x),a{x')) 



(1.4) 



The particular scaling guarantees that for the interaction 
potentials considered in the sequel the Hamiltonian re- 
mains finite as iV — > oo. To simplify notation we absorb 
the physical interaction parameters A, g, K into defi- 
nitions of corresponding interaction potentials. In the 
derivation of the statistical model we omit the external 
field to keep the notation simple and focused to the av- 
eraging procedure only. However, in the presence of an 
external magnetic field h,, the Hamiltonian also involves 
the interaction energy 



Uh 



(1.5) 



The central object describing the statistical ensemble 
on the lattice is the canonical Gibbs measure at the in- 
verse temperature P = -gf, defined on the configuration 
space S = {cr I a{x) e S*^, x £ /Cjv}, 



where 



Zn,p = / e-^^"('^)HAr(da) 



(1.6) 



(1.7) 



denotes the partition function. The prior distribution 
Hat (da) models the small scale fluctuations of the spins; 
we assume it is a product measure on the lattice Cj^ of 
uniform distributions on S^, i.e., that spins at different 
lattice sites are independent, uniformly distributed ran- 
dom variables on S'^ . Note that we can alternatively in- 
corporate the anisotropy energy in the prior distribution. 
Here we study the most probable configuration on Cn ac- 
cording to the Gibbs measure Pn,i3 ds N ^ oo. In this 
limit, we assume the lattice Cn approximates a physical 
domain H, C W^. We show that the energetically most fa- 
vorable configuration describing the large-scale features 
on Cn ~ is given by the Maxwellian distribution. 



M{x,v) 
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Zp{x) 



-l3[i!{v)+v{yu{x)-J*m(x)-h^)] 



(1.8) 



yielding the probability of having a spin u G S*^ at 
the location x G f2. Here Zp{x) denotes the partition 
function. The corresponding average magnetization is 
fh{x) — Jg2 vM(x,v) dv with magnetization potential u 
given by Au = div(xn?Ti). Furthermore (it, m) minimizes 
a newly derived free energy of Ginzburg-Landau type, at 
a finite temperature: 

Fg\m] = [ —iVul'^dx— [ —(J*m)mdx + 

a*p{m) dx — / h^rndx, (1-9) 



subject to the constraint 

Au = div(xnm) , in R'^ . 



(1.10) 



The equation ( 1.1C| ) , involving the characteristic function 
Xn of the domain Q, is understood in the weak sense (see 
Q). Furthermore J * m{x) = J{x — y)m{y)dy and 

a*p{m) = sup {mp - af3{p)} 
is the Legendre-Fenchel transform of the function 



(1.11) 



As in the case of the Boltzmann theory of dilute gases 
the equilibrium measure is defined by the macroscopic 
quantities: the magnetic potential u which defines the 
induced magnetic field h(x) — —Wu{x) and the mag- 
netization frL[x) — Jg2 vM{x,v) dv. A related approach 
is also adopted in the statistical description of coherent 
structures in 2D turbulence (see ||^). 

The equilibrium measure with the density M(a;, w) is 
identified as the minimizer of the energy functional 



Wul-'dx 



-(J * to) m,dx -|- 



n 



+ 11 ^{v) f{x,v)dvdx + 

+\ I I f{x,v)\ogf{x,v)dvdx. (1.12) 

Since the minimization is over the space of probability 
densities / the following constraints must be satisfied 



/(x, v)dvdx = 1 , 



together with ( 1.10 ) that relates u and m, where 



m{x) = I vf{x,v)dv. 

Ig2 



(1.13) 



(1.14) 



The last constraint in (1.13) is due to the fact that at 
every latt ice s ite there is one and only one spin. The 
quantity (1.14) evaluated at the minimizing density / — 
M, represents the average (macroscopic) magnetization. 
The last term of the energy functional is interpreted as 
the thermodynamic entropy and mathematically it is the 
relative entropy wit h re spect to the prior distribution of 
the Gibbs measure (1.6). 

We conclude the discussion of our main results with 
two brief comments on the random, thermally induced 
fluctuations. First, (1.8) is the probability density of the 



spin V at each spatial location x, and upon averaging over 
the spin space yields to. Furthermore, the large deviation 
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principle derived below with the rate Ep [f] , describes the 
equilibrium rando m fl uctuations around the most prob- 
able macro-state (1.8) as a function of the specimen size 
N. In the following sections we outline the main steps 
leading from the microscopic Hamiltonian to the macro- 
scopic variational principle involving energies Ep and Fa . 
The detailed mathematical derivation is described in [p| . 



II. THERMODYNAMIC LIMIT AND LARGE 
DEVIATION PRINCIPLE 

We employ the theory of large deviations (see [^|) to 
obtain the functional Ep from the spin Hamiltonian (1.4). 
We define the empirical measure on Af :— x corre- 
sponding to the spin configurations a, 



as N 



H [dy.dv) 



1 

TV 



x£Cn 



Integrating over the domain A x {i;} c M, we readily see 
that {A X {w}) is a coarse-grained random variable 
counting the number of spins v contained in the region 
A. In the N ^ oo limit it will give rise to the probability 
distribution of spins in A. The empirical measure allows 
us to show that, (i) the corresponding Gibbs states yield 
a coherent deterministic large scale picture, referred to as 
a "macrostate" and varying at the size of the domain 17 x 
5^, and (ii) the microscopic random fluctuations at the 
lattice cell size around the macrostates can be described 
explicitly in terms of Efj. 

A crucial step of the derivation is to compute a suitable 
limit of the lattice Hamiltonian and of the Gibbs measure 
(1.6) as N oo. Using the empirical measure we can 
write the lattice Hamiltonian as 



H{pi^) = - / U{x-x',v,v')n^{dx,dv)n^{dx',dv'). 

The partition function corresponding to the re-scaled 
Gibbs measure Pn,np is 



^N,Nfj 



S{p) 



Here 



f (x , v) log f (x , v) dxdv if dfi^fdxdv, 



and S{fj,) — oo otherwise, which is known in probability 
theory as the relative entropy with respect to the uni- 
form measure dxdv For the rigorous meaning of the 
approximation sign above in the framework of the the- 
ory of large deviations, we refer to ||^. By employing 
the asymptotics of = /x), as well as the Laplace 

principle, |9|, we compute the limit as oo 



lim ^ log Zn.n/s 



lim — loE 

N~>oo N 



-N{pH{p)+S(py) 



djjL 
(2.4) 



Similarly we obtain the large deviation principle for the 
Gibbs measure, as ^ oo, 

(2.5) 

The rigorous derivation is technically involved due to the 
nature of the singularity in the dipole-dipole interaction 
(see also ) . We address rigorously this issue in ||] . 

The expression ( |2.5| ) is interpreted as follows: the most 
probable configuration fi of the Gibbs measure is the min- 
imizeroi f3H{^) + S{p), yielding the large scale structure 
(at the size of the domain TV) at equilibrium. As we show 
below, this minimizer turns out to be the probability den- 
sity (1.8). Note here that /3iJ(/i) -I- S{^) is finite if and 
only if the measure /i has a density /, in which case we 
define the energy functional 

Ep[f]^H{f,) + ^S{t,). 

Substi tutin g for expressions on the right hand side and 
using (1.10 ) we immediately obtain (1.12). 

Finally, (2.5) captures the microscopic spatial random 
fiuctuations of the empirical measure /i^ as a function 
of N (i.e., the total number of spins in the specimen), 
around the mos t pro bable macrostate given by the min- 
imizer (O) of ([L12|). 



and using the Hamiltonian H we write 



1 

N 



logZ 



N,Nfj 



1 

N 



log 



where P is th e set of all probability measures jioxiM that 
satisfy ( 1.13| ) (observe that if /i^ converge in the sense of 



measures as iV — > oo to a probability meas ure o n J\f with 
density / then / satisfies the constraints ( I.13| )). Using 
Sanov's Theorem and an auxiliary coarse grained process 
(see [0) we can show the large deviation principle for 
the prior distribution with the rate Junctional S{fi): 



H^(M^=M)«e-^^(^') 



(2.3) 



III. VARIATIONAL PRINCIPLE 

We show that the Maxwellian M{x, v) which represents 
the equilibrium macrostate of the microscopic syst em i s 
the minimizer of the energy functional Ep. Using ( l.lOj ) 
we define an auxiliary functional. 



Wurdx 



-(J * to) mdx + 



+ 11 {^{v) +vVu) f{x,v)dvdx + 

15 f f f ix,v) log f{x,v)dxdv . (3.1) 
P Jn Js^ 
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We observe that E/slf] = sup^^ij^i^j^d-^ £p[f,u], and con- 
sequently 

inf [f] = inf sup £p [/, u] > sup inf £(3 [/, u] . 

f f u u f 

Note that the infimum is taken over the densities / that 



satisfy const raint s (1.13) and also relations (1.10) to- 



gether with ( 1.14 ). For a fixed u we can construct the 
minimizer of £p [f, u] over / by varying the probability 
density / using suitable push- forw ard maps (sec Jl2| , p^ ) 
that preserve the constraints ( 1.13 ). Then the minimizer 
M„, £^[M„, u] — minj iS/3[/, u] for given, fixed u satisfies 



Zp{x) 



'l3[^{v)+v(Vu(x)-J*'mu{x))] 



(3.2) 



where Zfi{x) is the corresponding partition function and 
fnu{x) = jg2vyiu{x,v) dv. Note that a^(Vu — J*TO„) = 
1/(3 log Zp {x) . Straightforward calculations show that 



muix) — dpapiJ * m„ — Vm) , 



(3.3) 



The existence of mu satisfying (3.3) follows from mini- 
mizing the functional 



apim) dx . 



1 f 

— {J*m)mdx+ / mS/udx 
n 2 Jq 



Hence we obtain that 

£j3[Mu,u] = min£p[f,u] = - ( -\Vu\^dx + 

+ / n {J * n^u) n^u dx — / apCVu ~ J * rriu) dx . 
Jn 2 Jo 

Using the duality between ap and a*^ we have 

sup inf £i3 [/, u] — sup £p [M„ , u] — iijf — . Fjs [ra] . 

u f u m&(1.10l 



Closer inspection of the functional Fp proves that the 
last infimum is attained and that a minimizer m satisfies 

m = dpafj{J * rfi — Vu) 



and m, u are related by (1.10). The expression ( |3.2| ) 
also implies the equilibrium Maxwellian M(a;,w) given 
by (1.8). Using the constructed fields rh{x), u(x) and 
M{x~vj we can check that, in fact, 

inf Ef3 [f] = inf sup 5^ [/, u] = sup inf £13 [/, u] = 

/ f U U f 

Fp[m] = Fp[fh] . 



m&(l.lC) 



Thus the minimizer of the new Ginzburg-Landau en- 
ergy i^/3[TO] provides the minimum value of the energy 
functional Ep[f] as well as it describes, through the 
Maxwellian density M(x, v) the structure of the most- 
probable macrostate. The rigorous treatment of all steps 
outlined in the previous sections is described in p]. 



IV. DISCUSSION 

A natural question is how the presented model relates 
to the Landau-Lifschitz theory ( jg,^ ) and how it can be 
interpreted in view of recent studies of domain formation 
based on Monte Carlo simulations ( ||l|,^). First we ad- 
dress the former issue; we expand the convolution J * m 
in the free energy ( |l.9| ), J * m — Jom + J2/2 Am + . . ., 
where Jq = fp J(r) dr, J2 — Ir]"^ J (r) dr , and after 



substituting into (l.£) we obtain 



Fp[m] - 



2 ' ' 



hwul'dx 



{a*p{m) - ^Imp) 



dx 



hi,mdx, (4.1) 



which can be interpreted as the finite-temperature ana- 
logue of the Landau-Lifschitz free energy: 



A. 



Vml dx 



•d L 



\IuV dx 



-f / '^(m)dx- \ he-rndx. (4.2) 
Jo Jo 

However, the free energy functional i^LL is minimized 
subject to the non-convex constraint |m| = 1. In the 
formulation presented here the direct calculation implies 
that since \dpap{p)\ < 1 we have a*p{m) — 00 for \m\ > 1 

and consequently the energy Fp (or Fp) is minimized 
subject to the constraint |m| < 1. This remarkable dif- 
ference from the Landau-Lifschitz theory is caused by the 
presence of the thermal fluctuations and the averaged na- 
ture of the magnetization m in our model. We conclude 
the discussion on the relation of our proposed finite tem- 
perature model to the Landau-Lifschitz theory with an 
interesting point: one would be tempted to say that a 
suitable zero-temperature limit (F-limit) of the free en- 
ergy Fp should yield the Landau-Lifschitz model or its 
relaxation. Indeed, the /3 — > 00 limit of the free energy 
Fp or Fp c an b e explicitly calculated, at least formally, 
yielding in ( |L9| ) and (^t|): 



a^(m) sup {/ : i?" i-^ i?, / < on 5^ ,/ convex} , 

if \m\ < 1 and a^(m) = 00, if \m\ > 1. We note that the 
same energy rel axat ion was obtained in [ p^ by a direct 
minimization of ( |4.2| ) over all admissible Young measures, 
when the energy exchange term ^/2|Vmp is neglected. 

Existence of magnetic domains and different types 
of magnetic walls in the Landau-Lifschitz theory is at- 
tributed to the competition of different contributions in 
the free energy and to the non-convex constraint |m| = 1. 
Due to the thermal agitation incorporated in the finite- 
temperature free energy Fp (or Fp) the norm |m| also 
fiuctuates. Although for low temperatures |m| can be 
close to the unit sphere, in general, we have the convex 
constraint |m| < 1. Nevertheless, our model allows for 
domain formation at finite-temperatures whenever the 
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exchange energy is strong enough compared to the tem- 
perature. 

More specifically, in the absence of an external field 
he, the condition that / — Jo9ppa^(0) is positive definite 
guarantees that m = and u = is the minimizer of Fp 
(as well as Fp) implying a uniform state, i.e., no domain 
formation. When the above condition is violated, then 
for a suitable domain Q and anisotropy, there exist non- 
trivial solutions m ^ and u = to the Lagrange-Euler 
equation for the energy Ffj (or Fp), 



(4.3) 



which have lower energy than the uniform state m = 0, 
M = 0. Similar steady states on 2D lattices, referred 
to as many-soliton solutions, were first analytically pre- 
dicted and also observed in Monte Carlo simulations in 
0. Under simplifying assumptions on the geometry of 
the domain it is not difficult to calculate the domain- 
wall profile, as well as its dependence on temperature, 
t he e xchange energy and anisotropy energy directly from 
(4.3). Note that since the magnetization m in our model 
is computed by averaging over the thermal fluctuations of 
spins the domains are identified with regions in O where 



m(a;)| « [ mg < 1 and rrip are non-trivial constant so- 
lutions of (4.3). In a forthcoming publication, domain 



formation and evolution, as well as the corresponding 
stochastic dynamics are studied both numerically and 
analytically in the context of the statistical theory de- 
veloped here. 
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